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Functionals on Triangulations of Delaunay Sets 



Nikolay P. Dolbilin' , Herbert Edelsbrunner*, Alexey Glazyrin^ and Oleg R. Musin* 



Abstract 

We study densities of functionals over uniformly bounded 
triangulations of a Delaunay set of vertices, and prove that 
the minimum is attained for the Delaunay triangulation if this 
is the case for finite sets. 

Keywords. Delaunay sets, triangulations, Delaunay triangulations, 
uniformly bounded triangulations, functionals, densities. 



1 Introduction 

A Delaunay set X C R d has positive numbers r < R such 
that every open ball of radius r contains at most one point, 
and every closed ball of radius R contains at least one point 
of X. Such sets were introduced as (r, /?)-systems by Boris 
N. Delaunay in 1924. By a triangulation of X, we mean a 
simplicial complex, T, whose vertex set is X and whose un- 
derlying space is R d . This triangulation is uniformly bounded 
if there is a real number q = q(T) such that the circumsphere 
of every ^-simplex in T has radius smaller than or equal to q. 
A particular triangulation is the Delaunay triangulation, de- 
noted as Del X, whose cZ-simplices satisfy the additional con- 
dition that all other vertices lie outside their circumspheres. 
It exists if X is generic, as will be explained shortly. The 
Delaunay triangulation of a Delaunay set is necessarily uni- 
formly bounded. We also consider Delaunay triangulations 
of finite sets of points, for which the underlying space of the 
simplicial complex is the convex hull of the points. 

Writing Sd for the set of all rf-simplices in R d , we con- 
sider functionals F : S c / — » R for which there are constants 
e = e(r, q, d) and E = E(r, q, d) such that e < F(cr) < E 
for every ^-simplex cr whose edges are longer than or equal 
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to 2r and whose circumsphere has a radius smaller than or 
equal to q. Writing & for this class of functionals, we de- 
fine subclasses Q c f c £ by requiring additional condi- 
tions. Briefly, F belongs to T if the sum of values over the 
(f-simplices of the Delaunay triangulation of d + 2 points in 
R d is smaller than or equal to the sum over the af-simplices 
in the other triangulation, and F belongs to Q is it satisfies a 
similar condition for all finite sets of points. For a triangula- 
tion, we define the density of the functional by taking sums 
and lower limits over a growing sequence of balls: 



f(T) = liminf 



1 



- Vol(B ) 



(1) 



where B a is the closed ball with radius a and center at the 
origin of W 1 , and Vol(B a ) is its volume. With these defini- 
tions, we can give our main result: 

• in R d , F e Q implies that the Delaunay triangula- 
tion minimizes the density of F among all uniformly 
bounded triangulations of a Delaunay set, and in K 2 , 
F € T suffices to reach the same conclusion. 

There are many concrete functionals studied in the literature 
to which our result applies. Here, we just mention two: 

• the functional that maps every triangle in R 2 to the ra- 
dius of its circumcircle; see O, 

• the functional that maps every rf-simplex to the sum of 
squares of its edge lengths times the volume; see IfTTI . 

The remainder of this paper presents the detailed results in 
two sections. 



2 Background 

In this section, we introduce the background on Delaunay 
sets, their uniformly bounded triangulations, and functionals 
on such triangulations. 



2.1 Delaunay Sets 

We recall from Section Q] that X c R d is a Delaunay set if 
there are positive constants r < R such that (I) every open 
ball of radius r contains at most one point of X, and (II) ev- 
ery closed ball of radius R contains at least one point of X. 
Hence, X has no tight cluster and leaves no large hole. 

Counting points. Condition (I) implies that every bounded 
subset of R d contains only finitely many points of X. In- 
deed, the subset can be covered by finitely many open balls 
of radius r, and each such ball contains at most one point. 
Condition (II) implies that every cone with non-zero volume 
contains infinitely many points of X. Indeed, the cone con- 
tains an infinite string of disjoint closed balls of radius R, and 
each such ball contains at least one point of X. We quantify 
the first observation by giving concrete estimates. Let B a (z) 
be the closed ball with radius a and center z, and call the 
difference between two concentric balls an annulus. 

Point Count Lemma. Let X be a Delaunay sets with pa- 
rameters r < R in R d , 

(i) There are constants p = p{R,d) and P — P(r,d) such 
that the number of points of X in B Qr (z) is between pa d 
and Pa d . 

(ii) There is a constant P' - P'(r, d) such that the number 
of points of X in B ff+ i(z) - B Qr (z) is at most P'a d ~ l . 

Proof. To prove (i), we note that B a (z) can be covered by 
some constant times (a/r) d balls of radius r, and that we can 
pack some other constant times (a/R) d balls of radius R in it. 
The lower and upper bounds follow. 

To prove (ii), we cover the annulus with a constant times 
a d-i i ^ b a u s f ra( ji us r. The upper bound follows. 

It should be clear that the bound in (ii) also holds for 
annuli of constant width, but not for annuli whose width is a 
positive fraction of the radius. 

Delaunay triangulations. Following the original idea of 
Boris N. Delaunay, we consider t/-simplices with vertices 
from X such that the open ball bounded by the (d- 1 )-dimension 
circumsphere contains no points of X. We call such <i-simplices 
empty. Here, it is convenient to assume that X is generic in 
the sense that no d + 2 points in X lie on a common (d - 1)- 
sphere. Under this assumption, the empty of-simplices fit to- 
gether without gap and overlap. Now consider two not nec- 
essarily empty but non-overlapping of-simplices that share a 
(d- l)-simplex, which is a face of both. Assuming the two d- 
simplices belong to a triangulation, we call this face locally 
Delaunay if the (d + l)-st vertex of the second (i-simplex lies 
outside the circumsphere of the first cZ-simplex. Note that the 
condition is symmetric because the two circumspheres inter- 
sect in the (d - 2)-sphere that passes through the vertices of 
the face, and either both (d + l)-st vertices lie outside or both 



lie inside the respective other circumsphere. Delaunay con- 
sidered both conditions and proved that they are equivalent 

m- 

Delaunay Triangulation Theorem. Let X be a generic De- 
launay set in R c '. 

(i) The collection of empty <i-simplices together with their 
faces form a triangulation of X, commonly known as 
the Delaunay triangulation, DelX. 

(ii) If all (d - l)-simplices of a triangulation T of X are 
locally Delaunay, then T = Del X. 

The equivalence between the local and the global conditions 
expressed in (ii) also holds for finite sets X. In the plane, it 
means that a triangulation of a generic set X is Delaunay iff 
for each edge the sum of opposite angles in the two incident 
triangles is less than n. 

2.2 Uniformly Bounded Triangulations 

Let X be a generic Delaunay set in R d , and let T be a tri- 
angulation of X. We recall that this means that T is a sim- 
plicial complex with vertex set X whose underlying space is 
R d . Recall also that T is uniformly bounded if there is a real 
number q = q(T) such that the radius of the circumsphere of 
every ^-simplex in T is smaller than or equal to q. It follows 
that no edge of T is longer than 2q. Note that the Delaunay 
triangulation of X is uniformly bounded with q — R. 

Not every triangulation is uniformly bounded. We begin 
by showing that every Delaunay set has triangulations that 
are not uniformly bounded. Given X, we construct such a 
triangulation in three steps. 

1 . For every point x e X and every L > 0, we can find 
many points y e X such that the edge xy is longer than 
L and does not pass through any other points of X. 
Indeed, there is such an edge near every direction out 
of x. To see this, we consider the set of points in X that 
lie within the closed ball of radius L around x. There 
are only finitely many such points, which implies that 
within each cone with non-zero volume and apex x, 
we can find a subcone, again with non-zero volume 
and apex x, that does not contain any of the points of 
X inside the ball. However, as argued in Section [2T1 
the cone contains infinitely many points of X, so they 
must all be at distance larger than L from x. Among 
these, let y be the point closest to x. 

Using the knowledge about long edges, we construct a tri- 
angulation inductively, one phase at a time. After the k-th 
phase, we will have a triangulation of a finite subset of X 
that includes all points at distance k or less from the origin. 
In addition, we will make sure that TV contains at least one 
edge longer than k, and that every point of X that belongs to 
the underlying space of is a vertex of T^. We start with 
To consisting of a single edge connecting the two points of 
X that are closest to the origin of R d . 
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2. In the (k+ l)-st phase, we let x be a vertex in the bound- 
ary of 7&. Let y be another point of X such that the 
edge xy is longer than k + 1 and does not intersect the 
simplices in TV other than at x. Since 7* is a triangula- 
tion, its underlying space is convex, and its boundary 
is triangulated. Let cr be an /-simplex in the boundary 
of 7* that is visible from y. We extend 7* by adding 
the (i + l)-simplex formed by y and the vertices of the 
/-simplex. Doing this for y and all visible simplices 
in the boundary of 7*, we obtain a simplicial complex 
T' k by starring from y. Similarly, we add a point z at 
distance k + 1 or less from the origin that lies outside 
the underlying space by starring to the triangulation. 
Repeating this operation for all such points z, we even- 
tually get a simplicial complex T'J. 

3. While 77' is a valid simplicial complex in E. d , some 
of the new simplices may contain points of X in their 
interiors. By construction, xy is not among these sim- 
plices. Let w e X be such a point, and cr e 77' the 
simplex of lowest dimension, j, that contains w in its 
interior. Note that j > 1. We fix the situation by de- 
composing cr into y + 1 y'-simplices, each the convex 
hull of w and j vertices of cr. Similarly, we decompose 
each simplex that contains cr as a face into j + 1 sim- 
plices of the same dimension by starring from w. Re- 
peating this procedure for all such points w, we even- 
tually get a triangulation Tk+\ such that all points in X 
that belong to the underlying space of T^+i are in fact 
vertices of T k+l . 

Observe that the edge xy added to 7* in Step 2 remains un- 
divided until the end. This implies that Tk+\ indeed contains 
an edge longer than k + 1, as required. It follows that the 
triangulation thus constructed by transfinite induction is not 
uniformly bounded. 

Measuring volume. The remainder of this section states 
and proves properties of uniformly bounded triangulations. 
We begin with the volume of their simplices. 

Volume Lemma. Let X be a Delaunay set with parameters 
r < R in R d , and let T be a uniformly bounded triangulation 
with parameter q of X. 

(i) In R 2 , there is a positive constant v = v(r, q) such that 
v < Area(cr) for every triangle cr in T. 

(ii) In R'', there is a constant V = V(q, d) such that Vol(cr) < 
V for every cZ-simplex cr in T. 

Proof. We prove (i) by expressing the area of a triangle in 
terms of the three edge lengths and the radius of the circum- 
circle: Area(cr) = The edges cannot be shorter than 
2r, and the radius cannot be larger than q, which implies 
v(r,q) = 2r 3 /q < Area( cr). 

To prove (ii), we note that every ^-simplex is contained 
in the ball bounded by its circumsphere. Since the radius is 
at most q, this ball is smaller than V(q, d) = (2q) d > Vol(cr). 



If we remove the requirement of uniform boundedness, 
then the proof of the Volume Lemma breaks down. It is not 
clear whether the upper bound fails. In this context, we men- 
tion a related question asked by L. Danzer and independently 
by M. Boshernitzen: "is it true that for every planar Delau- 
nay set there exists a triangle with arbitrarily large area that 
contains no points in its interior?" This question is still open. 

Next, we describe a Delaunay set in R 3 that has tetrahe- 
dra of arbitrarily small volume in the Delaunay triangulation. 
It shows that the limitation of the lower bound in (i) to two 
dimensions is necessary. Consider the standard cubic lattice, 
Z 3 . Let 5c = 2^j£|> f° r every t e Z, and move every point 
(i,j,k) e Z 3 to (i, j,k + (-1)' + ^), denoting the new point 
set by X. To study the volume of the tetrahedra in the De- 
launay triangulation, we consider a single integer cube, for 
which we get a tetrahedron of volume about | in the middle, 
four tetrahedra of volume about \ across each face, and two 

D 

flat tetrahedra at the top and the bottom; see Figure Q] These 




Figure 1 : A distorted cube decomposed into seven tetrahedra. The 
arrows indicating the distortion are exaggerated for better visibility. 

volume estimates assume arbitrarily small values of 6[ . If the 
third coordinates of the original vertices are k and k+1, then 
the volume of the top tetrahedron is l^+i, and that of the 
bottom tetrahedron is Since among the 5[ there are ar- 
bitrarily small numbers, there are tetrahedra in Del X whose 
volume is arbitrarily close to 0. It is not difficult to extend 
this example to four and higher dimensions. 

Counting simplices. Recall the Volume Lemma, which states 
that every triangle in a uniformly bounded triangulation of a 
Delaunay set in R 2 has an area that exceeds a positive con- 
stant. Since a disk of radius a has area a 2 n, this implies that 
the number of triangles contained in this disk is at most some 
constant times a 2 . A similar result holds in three and higher 
dimensions, but the lack of a lower bound on the volume of 
a rf-simplex requires a different argument, which we present 
as the proof of the following bounds. 

Simplex Count Lemma. Let X be a Delaunay set with pa- 
rameters r < R in E. d , and let T be a uniformly bounded 
triangulation with parameter q of X. 

(i) The number of <i-simplices sharing a common vertex 
is bounded from above by a constant S" = S "(r, q, d). 
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(ii) There are positive constants s = s(R, q, d) and S = 
S (r, q, d) such that the number of simplices contained 
in a ball of radius a > Aq is between sa d and Sa d . 

(iii) There is a constant S ' = S'(r,q,d) such that the num- 
ber of cZ-simplices contained in a ball of radius a + 1 
but not in the concentric ball of radius a is at most 

S'a d -\ 

Proof. To prove (i), we let x be the shared vertex, and we 
note that all incident li-simplices are contained in the ball 
of radius 2q centered at x. By the Point Count Lemma, 
the number of points in this ball is bounded from above by 
P(r, d) ■ (2q) d . We have at most one cZ-simplex for every com- 
bination of d of these points, which gives S"(r, q,d). 

The upper bound in (ii) is now easy: by the Point Count 
Lemma, the number of points inside the ball of radius a is 
at most P(r, d)-a d . Multiplying P(r, d) with S "(r, q, d) gives 
S (r, q, d). To get the lower bound, we restrict ourselves to the 
ball of radius a-2q. By the Point Count Lemma, the number 
of points in this smaller ball is at least p(R, q) ■ (a - 2q) d . 
Every ^-simplex incident to one of these points is contained 
in the ball of radius a. Each point belongs to at least d + 1 d- 
simplices, which implies that the lower bound on the number 
of points also applies to the li-simplices. Finally, (a - 2q) d is 
at least a d /2 d . 

To prove (iii), we use the upper bound of P'(r, d) ■ a d ~ l 
on the number of points in the annulus. Each of-simplex we 
count is incident to at least one of these points. Multiplying 
P'(r, d) with S "(r, q, d) gives S \r, q, d). 



2.3 Functionals 

Recall that Sd denotes the set of simplices, including degen- 
erate ones. We are interested in functionals that have con- 
stant upper and lower bounds for the simplices that arise 
in uniformly bounded triangulations of Delaunay sets. For 
other degenerate simplices we also allow infinity as a value. 

Definition. Let £ be the class of functionals F : Sd — > R 
for which there are constants e - e(r, q, d) and E - E(r, q, d) 
such that e < F(cr) < E for all af-simplices cr with edges of 
length at least 2r and radius of the circumsphere at most q. 

In this section, we extend the functionals from simplices to 
triangulations, and we introduce subclasses that favor Delau- 
nay triangulations for finite sets of points. 

Densities. As already mentioned in Section Q] we define 
the density of a functional on a triangulation by taking the 
lower limit over a growing ball, of the sum of values over all 
(f-simplices in the ball divided by the volume of the ball: 



f(T) 



lim inf „ , 

a^oo V0l(B o ) 



Z w - 



(2) 



There are other possibilities, such as taking the upper limit, 
or taking the average over the simplices. Our results extend 
to both modifications of the definition. One of Delaunay's 
motivations for defining (r, /?)-systems was to generalize lat- 
tices in Rr to a larger class of sets. For the Delaunay trian- 
gulation of any lattice A C W 1 , the limit of the expression 
in (f2]i, in which we substitute Del A for T, is well defined. 
Unfortunately, this is not generally the case for Delaunay tri- 
angulations of Delaunay sets, which is the reason for taking 
the lower limit. Since this might not be entirely obvious, 
we will prove shortly that for a broad class of functionals in 
£, the limit does not generally exist. Before that, we prove 
some positive results, namely that the density of every func- 
tional is bounded and independent of the choice of origin. 
Specifically, we define 



fz(T) 



lim inf ■ 



1 



Vol(B (z)) 



Z w 



(3) 



B„fe)2crE7" 



for every point z e R d , and we prove that all choices of z give 
the same result. 

Properties. Let F be a functional in £. 

(i) There is a constant C = C(r, q, d) such that f(T) < C 
for every uniformly bounded triangulation of a Delau- 
nay set in R d . 

(ii) f(T) = f : {T) for every z e R d . 

Proof. To prove (i), we recall the Simplex Count Lemma, 
which implies that the number of li-simplices contained in 
B a is bounded from above by S(r,q,d). Multiplying with 
E(r, q, d) gives C(r, q, d). 

To prove (ii), we let f(T, a) and f z (T, a) be the expres- 
sions in (f2j) and (0) without taking the lower limit, so that 
f(T) = lim info^co f(T, a), and similarly for f z (T) and f z (T, a). 
It suffices to prove 

lka[f(T,a)-f z (T,a)] = 0, (4) 

a— >oo 

which we do in two steps, namely by proving 

lim[f(T,a + L)-f(T, a )] = 0, (5) 

a— >oo 

]im[f(T,a + L)-f : (T,a)] = 0, (6) 

where L = \\z\\ is the distance between z and the origin. To 
prove (0, we write Vd for the volume of the unit ball in R d , 
and we note that Vol(B tJ ,) = Vda d ■ Furthermore, we write 
£fl and Ea for the sums over the tZ-simplices contained in the 
smaller ball and the extra <i-simplices contained in the larger 
baU: 



Z B = J] F(<r), 

B„2cte7" 



(7) 
(8) 
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By the Simplex Count Lemma, is at most some positive 
constant times a d , while is at most some constant times 
a d ~ l . Hence, 



A = f(T,a + L)-f(T,a) 



V d (a + L) d V d a d 

(a + L) d X B - a d X B 
V d (a + L) d V d (a + L) d a d 



(9) 
(10) 

(11) 



The first term in (fTTl i goes to zero because Sa grows slower 
than or, and the second term goes to zero because (a + L) d - 
a d grows slower than a d . The argument for © is similar. 
Indeed, all we need is to notice that the set of c/-simplices 
contained in B a+ L but not contained in B a (z) is a subset of 
those contained in B 0+ l but not contained in B ff _£. By the 
Simplex Count Lemma, the number of simplices thus defined 
is bounded from above by a constant times a d ~ l , so that the 
argument goes through as before. 



Non-existence of limits. We now show that taking the lower 
limit in the definition of density is necessary because the 
limit does not generally exist. Indeed, functionals for which 
the limit exists, even just for all Delaunay triangulations of 
Delaunay sets, are the exception. This is true in particular 
for the functionals that are invariant under isometries, which 
include all examples we discuss in this paper. 

We begin by exhibiting a construction in R 2 that acts as a 
stepping stone in our argument. Let cr and r be two compat- 
ible triangles, by which we mean that they share an edge, the 
two angles opposite that edge add up to less than n, and the 
remaining four angles are all acute. The condition implies 
that at least one of the triangles is acute, and we assume cr is. 
Using a linear sequence of congruent copies of cr, we form a 
strip T a , which we call wide, and using copies of t, we form 
a strip T T , which we call narrow; see Figure|2] Gluing strips 




Figure 2: A Delaunay triangulation made of a block of three wide 
strips in the center and two blocks alternating between narrow and 
wide strips glues on both sides. 

together so that they match up at boundary edges, we get a 
Delaunay triangulation, provided no two narrow strips are 
glued to each other. Let m\,m.2, ... be an infinite sequence 



of odd integers. We construct a Delaunay triangulation D in- 
ductively, starting with a block of tn\ wide strips. On each 
side, we add a block of mo strips alternating between narrow 
and wide, then a block of WI3 wide strips, then a block of tn.4 
strips again alternating between narrow and wide, and so on. 
For each i > 1, let 2a, be the total width of the first 2/ - 1 
blocks. Making sure that the origin lies on the center line of 
the first block, B a . is the largest disk centered at the origin 
that is still contained in the union of the first 2; - 1 blocks. 
We consider the sequence 



Si 



= — y 



F(o-). 



(12) 



B„.2o-£D 



Assuming that the limit in (f2]i exists, the sequence of f/ must 
converge for every sequence of m,. This is indeed the case if 
we measure area, because the number of triangles that inter- 
sect Bq. but are not contained in it is bounded from above by a 
constant times a. Hence, lim„^oo 4^ £ A(cr) = 1, where we 
abbreviate A(cr) = Area(cr) and take the sum over all trian- 
gles contained in B a , as usual. This motivates us to consider 
the ratios of the terms in the two sequences. Assuming B ffj 
contains k; congruent copies of cr and congruent copies of 
r, this gives 



DcreD 



A(cr) 



kjFjo-) + l,-F(r) 

M(cr) + £A(T) ' 



(13) 



Define go- = 77^, Q T = 77^, and Q = . Assuming 

Qa * Q T , we have Q a + Q and define A = \Q a - Q\. By 
definition, g\ = Q a . We choose m^ large enough so that 
l<?2 _ 21 < f > which is possible because ^2/^2 goes to 1 as 
1112 goes to infinity. Then we choose 1113 large enough so that 
|g3 - Qa\ < |, and so on, alternating between being close to 
Qcr and Q. We thus arrive at a contradiction because there is 
a gap of size | between the terms with odd and even indices. 



In other words, we need 



F(cr) _ F(t) 



for the limit to exist. 



A(tr) A(r) 

We finally show that the non-existence of the limit is 
not an artifact of the particular Delaunay set we used in the 
construction of D. Let cr' and t' be arbitrary triangles with 
longest edges of lengths a and c and opposite angles 2ip and 
2i/c Setting L > max L " , . c , }, we construct triangles cr 

" t> 1 2 cosy;' 2cost^' 

with edges of length L, L, a and t with edges of length L,L,c. 
It is easy to verify that cr' and cr are compatible, and so are 
cr and r, and r and t'. Repeating the construction with the 
strips and blocks three times, we see that if F is invariant 
under isometries of R 2 and the limit exists for the Delaunay 
triangulations of all Delaunay sets, then ^j^-j = Con- 
versely, among the functionals invariant under isometries, 
only the ones proportional to the area have the limit defined 
for the Delaunay triangulations of all Delaunay sets. 

Subclasses. We are interested in two subclasses of func- 
tionals, Q c f c £, which we now introduce. To define 
T, let Y be a generic set of d + 2 points in R rf such that 
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no point lies inside the convex hull of the others. The non- 
degenerate (f-simplices spanned by the points cover the con- 
vex hull twice; see Radon IflOl . Indeed, we can split them 
into two collections such that each forms a triangulation of 
Y: the Delaunay triangulation, D = Del Y, and the other tri- 
angulation, T. Changing one triangulation into the other is a 
flip, a name motivated by the planar case in which it replaces 
one diagonal of a convex quadrilateral with the other. We 
give the flip a direction, leading from T to D. Let now F be 
a functional, let be the sums of F(cr) over all of-simplices 
in T, and define Y,o similarly. 

Definition. The class T consists of all functionals F e £ 
for which So < £7-. 

In K 2 , the extra property of functionals in T suffices to prove 
our main result. In R d , for d > 3, we need more structure. 
The reason is the existence of triangulations that cannot be 
turned into the Delaunay triangulation by a sequence of di- 
rected flips; see for finite examples in R 3 . Such examples 
do not exist in R 2 ; see (7J. 

Let now Y be a finite set of points in TH d . As before, we 
assume that Y is generic. Let T be a simplicial complex with 
vertex set Y, but note that we do not require that T be a trian- 
gulation of Y . For example, we could start with a triangula- 
tion of Y and construct T as the subset of li-simplices that do 
not belong to the Delaunay triangulation together with their 
faces. Let D' be the subset of simplices in Del Y contained 
in the underlying space of T . Finally, let be the sum of 
F(cr) over all of-simplices in T' , and define T,jy similarly. 

Definition. The class Q consists of all functionals F e & 
for which £ D < < Y, T >. 

The condition for F to belong to Q is at least as strong as that 
for F to belong to T, which implies Q c f . 

3 Results 

In this section, we state and prove our Main Theorem and 
some of its implications. 

3.1 Main Theorem 

As already mentioned in Section Q] the main result of this 
paper is an extension of optimality results for Delaunay tri- 
angulations from finite sets to Delaunay sets, which are nec- 
essarily infinite. Section|2]provides all the technical concepts 
needed to give a formal proof of the theorem that facilitates 
this result. 

Main Theorem. Let X be a Delaunay set in R d . 

(i) In R 2 , F e T implies /(DelX) < f(T) for aU uni- 
formly bounded triangulations T of X. 

(ii) In W 1 , F e § implies /(DelX) < f(T) for all uni- 
formly bounded triangulations T of X. 



Proof. Fix a uniformly bounded triangulation T with param- 
eter q of X. We prove the inequalities by comparing subsets 
of rf-simplices of T and D = DelX. We begin with (ii). For 
every radius a, we write T(a) c T and D(a) c D for the 
sets of simplices contained in B a . Furthermore, we write 
D'(a) C D(a) for the set of simplices contained in the un- 
derlying space of T(a). Summing F over the cZ-simplices in 
these sets, we have 

^r(ff) = [?r(a) _ £/}'(«)] + [^D'(a) - ^D(a)] + ^D(a)- (14) 

The first difference on the right-hand side is non-negative 
by assumption of F e Q. We prove shortly that the sec- 
ond difference is bounded from above by a constant times 
a d ~ l . This implies that dividing by Vda d and taking the lower 
limit gives f(T) > f(D), as required. To prove the bound 
for the second term, we assume a > 4q, so that the Sim- 
plex Count Lemma implies that the number of of-simplices 
in D(a) - D(a - 2q) is bounded from above by a constant 
times a d ~ l . We get the same bound for D(a) - D'(a) because 
T(a) is uniformly bounded, with parameter q, so it covers all 
of B„_2fy. which implies D(a - 2q) c D'(a). 

The proof of (i) is similar, except that we have to do 
more work to construct the sets of fi-simplices, now trian- 
gles needed for the comparison. We assume a > 8g and let 
T(a) and D{a) be as before. We construct T'(a) by mod- 
ifying T(a) through a sequence of directed flips applied to 
non-locally Delaunay edges. Iterating the directed flip, we 
can guarantee that all interior edges of T'(a) are locally De- 
launay. A directed flip does not increase the size of the larger 
circumcircle (see e.g. 191), so flipping does not take us out- 
side the class of uniformly bounded triangulations. Impor- 
tantly, the flips turn a large portion of T(a) into Delaunay 
triangles, namely D(a - 6q) Q T'(a). To see this, we note 
that every triangle cr € T'(a) contained in B a ^ q is also in D. 
Indeed, its circumsphere is contained in fi«-2 9 , which is con- 
tained in the underlying space of T'(a). If cr were not empty, 
we would have a vertex inside the circumcircle, which would 
imply an edge that is not locally Delaunay between this ver- 
tex and cr, which is a contradiction; see also where this 
argument is used to prove part (ii) of the Delaunay Trian- 
gulation Theorem. Finally, the triangles of T'(a) contained 
in B a _4 9 cover B a -(, q , which implies D(a - 6q) c T'(a), as 
required. For the comparison, we consider 

^T(a) = Pr(a) - 2,T'(a)] + [^T'(a) ~ ^D(a-6q)] (15) 
+ [^Z3(ff-6<r/) ~ ^D(q-)] + ^D(a)- (16) 

The first difference is non-negative, and the second and third 
differences are bounded from above by a constant times a d ~ l . 
Dividing by V^ar 2 and taking the lower limit, as a goes to 
infinity, we get f(T) > f(D), as required. 

3.2 Implications in the Plane 

There are many functionals on triangles that are known to 
be in T ■ Applying the Main Theorem thus gives many opti- 
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mality results for Delaunay triangulations of Delaunay sets. 



Corollary A. Let cr be a triangle in R 2 , with edges of 
length a, b, c, let c\ > and C2 > 1 be constants, and consider 
the following list of functionals: 



• 


Fx{(f) = 


Circumradius Cl (cr). 


• 


F 2 ((r) = 


Circumradius C2 (cr) ■ Area(cr). 


• 


F 3 (cr) = 


-Inradius(cr). 


• 


F 4 (cr) = 


(a 2 + b 2 + c 2 )/Area(cr). 


• 


F 5 (cr) = 


(a 2 + b 2 + c 2 ) ■ Area(cr). 


• 


F 6 (cr) = 


||Centroid(cr) - Circumcenter(cr)|| 2 ■ Area(cr) 



Then fi(De\X) < fi(T) for every Delaunay set X c R 2 , for 
every uniformly bounded triangulation T of X, and for 1 < 
; < 6. 



Proof. It is easy to see that all listed functionals belong to £. 
For finite sets, the optimality of the Delaunay triangulation 
for f\ and f\ was proved in O, for fz and fe it was proved 
in J9), for fi, it was proved in J6), and for fi it was proved 
in HI1- It follows that F, 6 T for 1 < i < 6, so the claim 
follows from (i) in the Main Theorem. 



3.3 Implication in d Dimensions 

We have one example of a functional on <f-simplices that is 
in Q, namely the extension of F<, to three and higher dimen- 
sions. Writing a\ to ak for the lengths of the k = y^) edges 
of a cZ-simplex cr, we define Fr(ct) = Vol(cr) a 2 ; see also 
(TJ. Rajan proved that for finite sets in R d , the density of Fr 
attains its minimum for the Delaunay triangulation. We will 
extend his proof to show that Fr belongs to Q. With this, we 
get another consequence of the Main Theorem. 

Corollary B. We have /«(DelX) < /r(T) for every De- 
launay set X cR d and for every uniformly bounded triangu- 
lation T of X. 

Proof. The main tool in this proof is the lifting of a point y e 
R d to the point y + = (y, ||y|| 2 ) e R d+l , an idea that goes back 
to Voronoi 1121 . Note that y + lies on the graph of the function 
m : R d -> R defined by m(x) = \\x\\ 2 . For Y c R d , write Y + 
for the corresponding set of lifted points, and let conv Y + be 
its convex hull. Assuming Y is finite and generic, conv Y + is 
a convex polytope whose faces are simplices. We distinguish 
between lower faces whose outward normals point down - 
against the direction of the (d + l)-st coordinate axis - and 
upper faces whose outward normals point up. Importantly, 
if we project all lower faces vertically to R'', then we obtain 
the Delaunay triangulation of Y. 



Consider now the functional Fi that maps a of-simplex cr 
in R. d to the {d + l)-dimensional volume between the con- 
vex hull of the d + 1 lifted vertices and the graph of m. More 
precisely, it is the volume of the portion of the vertical (d+ 1)- 
dimensional prism over cr that is bounded above by the con- 
vex hull of the lifted vertices and below by the graph of m. It 
is not difficult to prove that Fi is invariant under isometries 
of R. d , and to use this fact to show that Fl € &. The reason 
for our interest in Fe is the relation 

F R (cr) = (d+l)(d + 2)F E (cr) (17) 

proved in IfTTI . Since the two functionals differ only by a 
multiplicative constant, it follows that Fr also belongs to S. 
It remains to prove that Fr belongs to Q, which we do by 
showing that Fe belongs to Q. Indeed, this should be clear 
from the lifting result: the lifted images of the of-simplices 
in the Delaunay triangulation are closer to the graph of rn 
than those of other J-simplices. More specifically, if T is a 
simplicial complex with finite vertex set Y in R d , and all sim- 
plices of D' c Del Y are contained in the underlying space of 
T , then the total (d + l)-dimensional volume we get for T 
is larger than or equal to that we get for D' . But this implies 
Fe 6 Q, and therefore Fr e @, as required. 



4 Discussion 

The main contribution of this paper is an extension of opti- 
mality results that hold for Delaunay triangulations of finite 
sets to Delaunay sets, which are necessarily infinite. In the 
plane, this extension holds for all functionals that improve 
upon flipping an edge that is not locally Delaunay. In three 
and higher dimensions, we need stronger properties to prove 
the extension. It would be interesting to know whether these 
stronger properties are necessary. Specifically, is it true that 
F € T implies that the density of F attains its minimum at 
the Delaunay triangulation of a Delaunay set in R d , also for 
d > 3? Similarly, are there functionals in Q that are not in 
f, or is T = @1 
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